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Abstract. 



1. Introduction. Topological degree [HIS] is one of the principal toolboxes of the 
modern theory of nonlinear dynamical systems. The range of applications of this 
toolbox is rapidly growing, see, for instance, [S]. In this paper we discuss a new, to 
the best of our knowledge, scheme of applying topological degree to the analysis of 
canard-type trajectories. 

If Vt^: K'' i-^- M'^ is a continuous mapping, J7 C M'' is a bounded open set, and 
y e R'' does not belong to the image W{dil) of the boundary dil of ft, then the 
symbol deg(T/F, ft, y) denotes the topological degree '4 of at y with respect to ft. 
If ^ W{dft), then the integer number j{W, ft) = deg(VF, ft, 0), called the rotation 
of the vector field W at dfi, is well defined. A detailed description of properties of 
the number "f{W, ft) can be found, for example, in P|. In particular, if id denotes the 
identity mapping, id(a;) = x, then the number 7(id — W, ft) measures the algebraic 
number of fixed points of the mapping W in ft. 
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We will use this tool to investigate canard-type periodic trajectories of singularly 
perturbed differential equations. Let us recall some related terminology. Consider 
the slow-fast system 



ez = h{x,y,z), 

where x, y, z are scalar functions of time, e is a small positive parameter, and /, 
(7, h are scalar functions. The subset 



of the phase space is called a slow surface of the system ([T]): on this surface the 
derivative i of the fast variable is zero. Moreover, a part of S where 



is called attractive {repulsive, respectively). A line L <Z S which separates attractive 
and repulsive parts of S will be called a turning line. In what follows, we suppose 
the turning line to be smooth. Trajectories which at first pass along, and close to, 
an attractive part of S and then continue for a while along the repulsive part of S 
are called canards or duck-trajectories [2 . 

In the paper, we focus our attention on periodic canards in singularly perturbed 
picewise linear systems, or more specifically, in a special case of ^ with h[x, y, z) = 
x-'r\z\. This choice has a number of reasons. From the methodological point of view, 
piecewise linear systems are convenient because they are integrable. Furthermore, 
such systems are used extensively in the modelling of a wide range of physical 
processes and have applications in electrical circuits [71 [31 E], flight control [H], 
chemical processes control [13] and neural subsystems with control behavior [S]. 
The piecewise linearity may be due to nonlinear elements such as saturation or 
may result from linearization about various operating points of a nonlinear plant. 
For example, the McKean model is a piecewise linear caricature of the FitzHugh- 
Nagumo model [10] . 

The main mathematical reason to consider the existence of periodic canards for 
this class of differential systems is that they (together with classical canards) are 
mathematical explanations of the limit behaviour of periodic solutions of a family of 
differential equations. However, the traditional methods of the "chasse au canard" 
are adapted for sufficiently smooth systems, and we consider the use of topological 
methods in the case of Lipschitzian nonlinearities as a dire necessity (continuous 
piecewise linear functions satisfy the Lipschitz condition, see, for example, [6]). 

Note that system ([T]) has a vector slow variable {x,y), and it is known that the 
existence conditions for canards in systems with scalar and vector slow variables 
have a vital difference. In the scalar case it is necessary to have an additional 
parameter in the system under consideration, which is demonstrated in [151 112] 
where piecewise linear systems on the plane are studied. The canard then exists 
only for a small range of values of the parameter. In the vector case such a parameter 
is not required, hence the system ([T]) , under some assumptions about the right hand 
side given in the next section, always has a canard. Below we show that ^ also 
has a periodic canard if some additional assumptions are made. 



X 



y 



f{x,y,z), 
9(x,y,z), 



(1) 



S* = {(x, y, 2) e : h{x, y, z) = O} 
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2. Main result. In this section we formulate the main existence result for topo- 
logically stable canard type periodic trajectories, using a simple example. 
Consider a system of ordinary differential equations: 

i = f{x,y,z), 

y^9(x,y,z), (2) 

£Z = X + \z\, 

with the additional assumptions described below. 

The slow surface of the system consists of two half-planes, an attractive half- 
plane Pa and a repulsive half-plane Pr : 

Pa = {{x,y,x):x<0}, (3) 

Pr^{{x,y,-x):x<0}, (4) 

together with the turning line 

L = {(0,y,0)}. (5) 
We consider auxiliary equations 

i ^ fa{x,y) ^ f{x,y,x), 
y = 9a{x,y) = g{x,y,x), 

and 

i = fr{x,y) ^ f{x,y,-x), 
y = 9r{x,y) = g{x,y,-x), 
which describe the dynamics near the slow half-planes ([3|) and (jj]) in the limit e — > 0. 

Assumption 1. The functions f and g in the right-hand side of ([2]) are globally 
bounded and globally Lipschitz continuous with a Lipschitz constant X: 

\f{x,y,z)l\g{x,y,z)\< M. (8) 

Assumption 2. The following relationships hold: 

/(0,0,0) = 0, 3(0,0,0) >0, y/(0,y,0)<0. 

An important role is played below by the solutions of ([6]) and ([7]). Denote by 
^a(*) = iXa{t),yait)) the solutiou of the system ([6]), satisfying the initial condition 
x{0) = 2/(0) = 0, and by w*{t) = {x*,y*) the solution of the system d?]), satisfying 
the same initial condition. 

Assumption [5] ensures the existence of a finite time interval {Ta,Tr) 3 such 
that 

xlit) < for Ta < t < 0, x;{t) < ioT < t < Tr. 

Assumption [2] also implies strict limitation on the possible location of canards of 
system ([2]), which should follow closely the attractive half-plane ([3]) for a certain 
interval ta < t < 0, and then move along the repulsive half-plane ([4]) for < i < ^r- 
Any such canard must first follow the curve 

^a = {{x:{t),y:{t),x:mcPa (9) 

for negative times, and then the curve 

Tr = {{x:{t),y:{t),-x:{t))} c Pr (lo) 
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for positive times, passing near the origin, where the two curves meet att — 0. Using 
standard tools, see [HITl], it is easy to see that for any time interval [ta, t,.] C {Ta,Tr) 
with ta < < tr such a canard exists for every sufficiently small e > 0. 

The question whether there exist periodic canards is less obvious. The above 
argument shows that a periodic canard should have a segment of fast motion from 
a small neighborhood of some point of the curve to a small neighborhood of the 
curve Ta] this fast motion is, consequently, almost vertical (i.e., almost parallel to 
the z axis). More precisely, if there is a limit of periodic canards as e — 0, then 
the limiting closed curve has necessarily a vertical segment connecting and F^. 
The next assumption ensures a possibility of such vertical jumps. 

Assumption 3. The trajectories w*(t) andw*(t) of systems ([6]) and ([7]) intersect, 
that is, there exist r and a such that 

x*a{T) ^ x*{c7) ^ X* , y*a{T) = y*{a) = y* 

with 

Ta <T <0 < a <Tr. 
Assumption 4. The intersection is transversal, that is 

A = fa{x*,y*)grix*,y*) - fr{x\y*)ga{x\y*) ^ 0. 
Assumptions [3]-[4] are illustrated on Fig. [TJ 




Figure 1. Solutions and w*{t). 

The following theorem states the existence of periodic canard in the system ([2]) 
under the Assumptions [T]-|4l 

Theorem 1. For every sufficiently small £ > there exists a periodic solution of 
system The minimal period Tniin of this solution approaches a ~- t as e — >■ 0. 

It will also be shown that the periodic canard passes through a small neighbour- 
hood of the point (cc*, j/*, 0), with the diameter of this neighbourhood going to zero 
as e goes to zero. 

Below we discuss the main ideas behind the proof of this theorem. The proof 
itself, which is divided into several technical propositions and lemmas, will be pro- 
vided in the next section. 
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Denote by 

Wa{t,to,xo,yn) = {xa{t,to,xo,yo),ya{t,to,xo,yo)) 
the solution of system ([6]) with the initial condition 

x{to) = xo, y{to)^yo- 

Similarly, denote by 

Wr{t,tQ,xo,yo) = {xr(t,to,xo,yo),yr{t,to,xo,yQ)) 

the solution of system ^ with the same initial condition 

x{tQ) = xo, y{to)^yo. 

Consider a small vicinity of the point {x* , y*). Since the intersection between Ta 
and is transversal according to Assumption SI the numbers u{x,y) and v{x,y) 
may be defined in this vicinity by 

Wa{u{x, y),0, X, y) G T^, Wr{v{x, y),0, x, y) G Tq. 




Figure 2. The set U{a) on the plane {x,y). 

Using the coordinates u(x,y) and v{x,y) we introduce for a sufficiently small 
a > a 'parallelogram' set n(a), illustrated on Fig. [2| 

n(a) ^{x,y: \u{x,y)l\v{x,y)\ < a/2}. (11) 

Also introduce the notation for the two 'sides' of this parallelogram: 

Q- = {x,y- v{x,y) = +Q:/2sgnyl, \u{x,y)\ < a/2}, 

Q+ = {x,y- v{x,y) = -a/2sgnA, \u{x,y)\ < a/2}. 

Denote by 

We{t,XQ,yQ) = {xe{t,xo,yo),y^{t,xo,yo),z^{t,xo,yo)) (12) 
the solution of the system ([2]) with the initial condition 

x{t) = xo, yij) = 2/0, z{t) = 0, 

with (a;o,2/o) G n(a). 
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In our assumptions, the solution We{t, xo,yo) first rapidly approaches the attrac- 
tive half-plane Pa defined in ([3]), and then follows Pa until t 0. The possible 
subsequent behaviors of the solution are the following: 

1. The solution stays close to the attractive part of the slow surface for t < 6, 
where 6 > 0. 

2. The solution begins to rapidly increase in the z direction around t « 0. 

3. The solution follows the repulsive half-plane Pr defined in ([4]), until a moment 
t* . Moreover, for positive t < t* it follows the curve F^. After the moment 
t* the solution may begin to rapidly increase in the z direction, or it may fall 
down back to the attractive part of the slow surface. 

To distinguish between these cases, we define the time moment 

Se{xo,yo) = min{{Tr} U {t: i{xo,yo) + p <t <Tr: Ze{t,XQ,yQ) < 0}), (13) 

where t{xo, yg) is a time moment close to zero, which is associated with the solution 
T, xq, 2/o) E^iid which will be introduced in the next section, and p > is a 
sufficiently small number chosen together with a. 

Using the moment s^ixo, yo), we define an auxiliary mapping of the parallel- 
ogram n(a) into the plane {x,y). The definition is divided into the following four 
cases: 
Case 1. If 

a ~ a < Se{xo,yo) < cr + a, 

then 

We{xo,yo) = ixe{se,xo,yo),ys{ss,Xo,yo)). 
If we identify the plane {x,y) with the two-dimensional subspace 

Po^{ix,y,0):x,yeR} 

of the phase space of system ([2]), then in Case 1 the value We{xo,yo) coincides 
with the intersection of the trajectory (fT2|) with Pq, as long as the corresponding 
intersection time is close to cr. 
Case 2. If 

a < \se{xo,yo) - a\ < 2a, 

then 

Tir / \ 2q; — Isg — crl , , , , Is — al—a^,. 
We{Xo,yo) = {Xe{Se,Xo,yo),ye{SE,Xo,yo)) H W^(Ss). 

a a 
This means that Wf;{x(), j/o) is a continuous convex combination of the intersection of 
the trajectory (jl2p with Pq and the point w*{si;) as long as the discrepancy js^ — cr| 
between the corresponding intersection moment and a in the the range from 
a to 2a. Moreover, if the discrepancy equals a, then the definition is consistent 
with Case 1; if the intersection moment equals cr ± 2a, then coincides with 
w;{a±2a). 

Case 3. If 

Se{xo,yo) >(J + 2a, 

then 

WeixQ,yQ) = w*{a + 2a). 

Case 4. If 

Ss{xo,yo) < cr - 2a, 
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then 

We{xQ,yQ) = w*{a - 2a). 

For small e, a and p the time moment Se{xo, yo) depends continuously on (a;o, J/o), 
therefore the operator We(xo,2/o) is also continuous with respect to (xo^yo). We 
will show that if (a;o,yo) is a fixed point of then Case 1 takes place, so the 
triple (xq, 2/0,0) defines a periodic solution of ^ with the period — r close to 
a — T. 

We also prove that the images of the parallelogram sides (3_ and (5+ are the 
points ■w*{a + 2a) and z/;*(cr — 2a) correspondingly, and u{We{xo,yo)) — J> as £ — 0, 
thus the operator is contracting along the u coordinate, and stretching along 
the V coordinate. Finally, the vector field rotation theory is used to obtain the 
relation 

7(id- We,n(a)) =sgnA^O, 

which implies the existence of a fixed point of the operator We in the set n(a), 
thus proving the existence of a periodic canard. 

3. Proof of Theorem [Tl The proof consists of several parts. First, we formulate 
the properties of the solutions Wa(i,xo,yo) = Wa{t,T,xo,yo) of the system ([6|) on 
the set Il{a). Secondly, we consider the behavior of the solution We{t, xq, yo) of the 
system Q and establish the continuity of the time moment Se{xQ,yQ) which was 
defined in (I13p and thus the continuity of the operator on n(Q;). Finally, we 
calculate the rotation of the vector field id — on 11 (a). 
We will need the following auxiliary definitions. 

Definition 1. We say that a is e-close to 6, if a — 6 as £ — > 0. a is said to be 
s-small, if it is e-close to 0. 

Consider a solution Wa{t, XQ^yo) of the auxiliary system ([6]) with the initial con- 
dition {x{T),y{T)) = {xo,yo) G 11. 
Define the set 

R = {{x, y): {x = OAy<0)V{x<OAy^ 0)}, 

which is the union of the left half of the horizontal x coordinate axis and the bottom 
half of the vertical y coordinate axis. Consider all the time moments T{xo, yo) when 
the solution Wa{t,Xo,yo) intersects the set R: 

T{xo,yo) = {t: Wa{t,xo,yo) G R}. 




Figure 3. The set T and the moment i. 
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Definition 2. Let i{xo,yo) be the moment from the set T{xo,yo) closest to zero, 
see Fig. [3] 

t = argmin \t\. 

tGT(xo,yo) 

Below we sometimes omit the point (xo,yo) and write simply t, in which case the 
arguments can be uniquely identified from the context. 

Definition 3. If t{xo,yo) exists, and Xa{t, xo,yo) — and ya{t,xo,yo) < 0, then 
the solution Wa{t^xo,yo) together with the point (xo,2/o) are called destabilizing. If 
Xa(t, xo,yo) < and ya{t,xo,yo) = 0, then the solution Wa{t,xo,yo) and the point 
{xotUo) are called stabilizing. 

This definition emphasizes the fact that the solution Wi;(t,xo,yo) of the main 
equation ([2]) with a stabilizing initial condition will stay close to the attractive 
half-plane Pa for some time after i > 0, if £ is sufficiently small; if initial condition 
is destabilizing, then Zi;{t) will rapidly increase to infinity after t > 0. The proof of 
this fact will be the subject of several propositions, all leading to Proposition [HI 

Statement 1. If a is sufficiently small, then the moment t is defined for all 
{xotUo) £ n and is continuous on 11 with respect to (xQ^yo). 

Statement 2. If a is sufficiently small, then a point (xo,?/o) £ 11 is destabilizing if 
A-v{xo,yQ) > 0, and stabilizing if A-v{xQ,yo) < 0. In particular, is destabilizing 
and (5_|_ is stabilizing. 

Statement 3. Let a be sufficiently small. If{xo,yo) destabilizing, thenXa{t) < 
for T < t < t, and if (xq, j/o) stabilizing, then Xa{t) < for t < t < (5(a) where 
(5(a) > maxn{t(a;o, yo)} > depends only on a. 

Statement [1] follows from the continuous dependence of Wa on {xo,yo) and from 
the fact that w* intersects the line y = transversally. 

To prove Statements [2] and |3l we can consider the projection d{t) of the difference 
between the trajectories Wa{t, xo,yo) and w*{t) onto the normal vector for w*{t), 
t > t: 

d{t) = {xa{t) - x:{t))ga{x:it), ylit)) - {ya{t) - y^t)) f^ix^t) , yHt j). 

Then the variation of d{t) satisfies the following initial value problem: 

m = {-§-Ja{xl{t),y:{t)) + f^9a{xl{t),y:{t)))r{t), r{T) = 1, 

and the following equality holds for small Axq = xq — x^,, Aj/o — yo — y*, = 
Axoga{x^,y^) - Ayofa{x^,y^): 

d{t) ^ r{t)do + o{Axo) + oiAyo). 

Therefore, if do > 0, and Axq and Ayo are small, then d{t) > for r < t < T^. 
Assumption!?] provides that do > for {xo,yo) G Q-, and do < for {xo, yo) G Q+- 

Denote a; = e In — with T = —Ta + T^., w > e for sufficiently small e. 

Consider the solution Wi;{t,xo,yo) defined by (jl2p with {xo,yo) £ n. 

Proposition 1. The sets z > —x + Me and z > x — Me are invariant in the sense 
that once a trajectory enters them, it does not leave them. In particular, if xq < 
and zo — 0, then > — Me for all t > t, and if additionally ^^(to) > —x^ + Me 
for some to, then z^(t) > for all t > to. 
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Proof. Denote ip{t) = z{t) + x{t), ■0(0 = z{t) - x{t). Then ^ together with 
assumption ([8]) imply that 

= i + i = -ix + \z\) + f(x, y, z) > -{x + z)~ M ^-(p-M, 
e e e 

tP = z-x = -{x + \z\) - f(x, y, z) > -ix - z)- M = --tfj - M, 
e e e 

and by the theorem on differential inequalities, 

ip{t) > ^{t) = (^0 - Me)e^ + Me, 
^{t) > 4'{t) = (V-o + Me)e~*-^ - Me, 

with LpQ — (y5(<o) — z{to) + x{to). and V'o = V'(^o) — z{to) — x{to). If ipo > Me, then 
ip{t) > Me and thus z{t) > —x{t) + Me for all t > to. Similarly, if tpo > —Me, then 
ip{t) > -Me, therefore z{t) > x{t) - Me for all t > to- 

The second part of the proposition follows from the fact that the point (xq, ya, 0) 
with xo < is in the set z > x — Me, and if both z > x — Me and z > —x + Me 
hold, then z > 0. □ 

Proposition 2. Let e be sufficiently small. Denote 

i = min({T,} U{t>T: z,{t) = 0}), (14) 

and 

,,.. + .l„(-i^L)^ ,15) 

Then ti < t + uj, t > ti, and for ti < t < t: 

\zeit) - x,{t)\ < 2Me. 

Proof. According to Proposition[TJ the solution w^{t) does not leave the set z — x > 
—Me, because {xq, yo) € 11 and thus zq- xq > 0. Therefore, it suffices to prove the 
relation Zs{t) - Xe{t) < 2Me for ti <t<i. 

First we prove that z^(t) < for r < t < ti, which will imply t > ti. Note that 
—MT < xo < —Me, thus ti < t + w. Assuming that z{t) < 0, we get the following 
relations for ip = z — x: 

^ = z-x = -{x- z)- f{x, y, z), 

e 

ij > - M, ipo^ zo-xo> Me, (16) 

e 

^(t) > {iPo + Me)e~i'^*-''^ - Me = V'^HO- 
The function ■ip'-^^{t) is decreasing, and 

4'^^Hti) = (^Ao + Me)— - Me > 0, 

therefore = —-tp < —-^p'-^'> < for r < i < ii, which proves that z^{t) < for 

T<t<tl. 
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Now, similarly to (jl6p . we write 

ip < --ip + M, 

e 

^{t) < {^0 " Me)e-^^(*-^) + Me = ^j'^'^Ht), 
^^^\tx) = (^0 - Me)^ + Me < 2Me, 

thus for t > ti: 

zAt)-Xe{t) < V'^^n^) < 2Me. 
This inequality holds while z^{t) satisfies the equation z = i(x — z), i.e. while 
t<i. □ 

As before, {xa,ya) will denote the solution of the system ([6]) with the initial 
condition x{t) — xq, y^r) — yo- 

Proposition 3. // e is sufficiently small, then the following inequalities hold for 
T <t< i: 

\Xe{t)-Xait)l\ye{t)-ya{t)\ < LlV , 

where L = bMe^^'^ . 

Proof. The moment ti defined by p5|) in Proposition [2] is w-close to r, thus 

\xeiti) - Xa{ti)\ < \xeiti) - xo\ + \xaiti) ~ xa\ < 2Muj, \y,{ti) - yaiti)\ < 2Muj. 

Denote il){t) = \xe{t) - Xa{t)\ + \ye{t) - ya{t)\, t > ti. Using ^ and © together 
with Assumption [T] and Proposition [21 we obtain 

DfilXe - Xa\ < |/(a;e,ye,Xe) - f{Xa,ya,Xa)\ + 2XMe < 2X\Xe - Xa\ + X\ye - J/a | + 2XMe, 
Dnlye -ya\ < \gixe,ye,Xe) - g{xa,ya,Xa)\ +2XMe < 2X\xe -Xa\ + X\ye -ya\ +2XMe, 
where Dn denotes the right derivative. Therefore 

DB.^j{t) < 4AVXt) + 4AM£, ij{ti) < AMoj. 
and by the theorem on differential inequalities, 

iPit) < (4Mcj + M£)e''^(*"*i) - Me < Luj. 
for ti <t<i. □ 

Proposition 4. // (a;o,?/o) destabilizing and e is sufficiently small, then the 
solution {x^,ye, z^) reaches z = at the moment t which is e-close to t{xQ,yQ), and 
Zeit) > for all t > i. 

Proof. Denote y = y{i) < 0, f — /(0,?/,0) > 0, .g = g{0,y,0). Consider a small 
vicinity ft of the point (0, y, 0): 

^ = {{x,y,z): \x\, \z\, \y-y\< S} 

Due to the continuity of / and g we can select a sufficiently small S such that 
f/2 < f{x,y,z) < 2f, and \g{x,y,z)\ < 2max{/, Lg]} for ah {x,y,z) G n. Denote 

At = V(2niax{/, I5I}), i - .x„(t- At) < 0. 

Then {xa{t),ya{t),Xa{t)) G for t — At < t < t + At, Xa{t) is increasing on this 
interval, and Statement [3] implies that 

a;Q(t) < i for T < t < t - At. (17) 
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First we will prove that t is e-close to i. Let At — ALui/f < At for sufficiently 
small e. We have 

Xaii- At) < -2Luj, Xa{t + At) > 2Luj. (18) 

Consider the interval I = {t: ti <t< t — At}. Inequalities ([T7|) . ([TSl) together with 
the fact that Xa{t) is increasing for i — At < t < i — At imply that Xa{t) < —2Luj 
for t ^ I. Applying Propositions [2] and O we obtain 

z^{t) < Xeit) + 2Me < Xa{t) + Luj + 2Me < for i e /, 

therefore z^{t) <{) ior t <t <i ~ At. 

Now suppose that Ze{t) <Q ioT i — At <t <i + At. Then Propositions [2] and [3] 
will hold on this interval, and (jl8|) will imply that 

ze{i + At) > Xe(t + At) - 2Me > Xa{i + At) - Llu - 2Me > 0. 

This contradiction proves that i lies in the interval i — At<i<i + At. 
Finally, we need to prove that z^{t) > for t > i. We have 

z,{i) = -{x,{i) + z,{i))>o, 

e 

thus Xe{t) > 0. It will suffice to show that the solution enters the set z > —x + Me 
before the time i + 0{e) (so that the solution will still be in the set il by that 
time, and Xi;{t) will be increasing, which will guarantee that z^{t) > 0). Let ip(t) — 
Zeit) + Xs{t), then 

^(t) = -^it) + f{x,,y„z,) > -^(t) + i/, ^{i) > 0, 
e e z 

and by the theorem on differential inequalities. 



<y5(0 > ^fe (e*(*-*) - l] for t > t. 



Consider the time moment i + e\n . At this moment 

4Af\ 1 ^ [AM \ 1 
— > _ f p — ' 



The proposition is proved. □ 

Proposition 5. // (xq, yo) ^•s stabilizing and e is sufficiently small, then < for 
all T < t < (5(a). 

Proof. According to Statement [3l Xa{t) < for t < i < 5{a). Denote 

X = min Xa(t) < 0. 

T<t<S{a) 

If e is sufficiently small, then Xa{t) < ~2Me on this interval, thus 2;e(t) < Xe{t) + 
2Me < for aU ti<t< S{a). □ 

Denote 

Se(xo,yo) = min({rr} U {t > i{xo,yo) + p: Ze{t) < 0}). 

The parameter p is sufficiently small and is chosen in the following way: due to 
Assumption [21 there exists a (5- vicinity of (0,0,0) such that f{x,y,z) is small 
and g{x,y,z) is positive. Then p must be less than 6/{2M), so that if a solu- 
tion (xe,?/e,Ze) is e-small at the moment i, then it remains in this vicinity for 
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i — 2p<t<t + 2p, and increases on this interval. Additionally, p should satisfy 
the inequality p < a/4. 

Proposition 6. // [xq, yo) is destabilizing, then for sufficiently small e, s — Tr- If 
(xq, j/o) stabilizing, then for sufficiently small e, s — p. 

Proof. Follows from Propositions |4] and [5) □ 

Proposition 7. Let Ze{s) > and e be sufficiently small, and let t be given by 
(|T4l) . Then \i-i{xo,yo)\ < p/4, < \xe{i)\ < Me, and ye{i) is e-small. 

Proof. From the definition of t it follows that z^{t) — ■i(a;e(t) + > 0, thus 

Xe{t) > 0, and the fact that z^{s) = implies that the solution never enters the set 
z + X > Me, which is invariant according to Proposition [1] Therefore Xe{t) < Me. 

Denote ye{t) = ij- Using the same reasoning as in Proposition 21 it can be shown 
that whatever y < is, for sufficiently small values of e the solution z^{t) > 
for t > t. Similarly, if y > 0, then f{x,y,z) < in some vicinity of the point 
(0,y, 0), and for sufficiently small e there exists At such that a;e(t — At) > 2Me, 
z^{t — At) < 0, which contradicts Proposition [21 

Therefore, there exists a function 7(e) as e — 0, such that \y\ < 7(e). 

It remains to prove that |t — t{xo, yo)\ < p/A. First we show that t > t{xo, yo) ^ 
p/A. This follows from the fact that Xa{t) < for t < t < t — p/4, so for sufficiently 
small e, x^it) < Xait) + Le < for t < t - p/4. 

To prove the relation t < i + p/4, we note that Proposition [51 implies 

\xaii)\,\yaii)\<i{e). 

Suppose that t > t + p/4. Then, taking into account Propositions [H and [31 we get 
that both j/e(t) and ye{t) are e-small, which is impossible because in this case y^ 
should increase for t < t < t due to the choice of p. Therefore, t < t + p/4 and the 
Proposition is proved. □ 

Proposition 8. Let z^{s) > with a sufficiently small e. Denote 

t = min{t < s: z^{^) >Ofort<^<s}. 

Then |t — t{xo,yo)\ < p/2, < |a;e(t)| < Me, and ye{t) is e-small. 

Proof. By the definition of s and t, t < t < i + p. According to Proposition [71 
the solution is e-small at the moment t. Thus, due to the choice of p, the function 
ye{t) increases on the interval t < t < t. Therefore, ye(t) > ys{t) > —7(e) with an 
e-small 7(e). 

Consider the solution {x^,y^, z^) in backward time at the moment t. Using the 
same technique as in Propositions |4l and |7l it can be shown that there exists an 
e-small function 7(e) such that yeit) < 7(e)- 

It remains to show that i < i + p/2. This follows from the fact that t < t -I- p/4, 
thus the function ye increases for t < t < t , therefore the e-smallness of both ye{t) 
and ye(t) implies that t and t are also e-close. □ 

Proposition 9. Let Zg(s) > with a sufficiently small e. Then 

\ze{t)+Xeit)\<2Me (19) 

for t < t < s — uj, and 

\xS) ~x:it-t)\, \yS) - y*r{t -t)\< i7, (20) 
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for t < t < s, where 7 is s-small. 

Proof. To prove the formula (jl9p . consider the solution {xe,ye, Ze) in backward 
time with the initial condition {xe{s),ys{s), Zs{s)) at the moment s. The the proof 
follows that of the of Proposition [21 with the exception that the fact that Ze{t) > 
for s — ijj<t<s is now provided by the statement of the Proposition {s — i > p/2, 
thus the interval above is non-empty). 

The proof of the formula ([^0)1 uses and Proposition [51 and is similar to that 
of Proposition [3l □ 

Proposition 10. If e is sufficiently small, then the moment of time Se{xo,yo) is 
continuous on H with respect to xq and j/o- 

Proof. Let {xo,yo) G n and prove that s is continuous at the point (xo,?/o)- There 
are three possible cases: 

1. Zf;{s) < 0, s — t+p. Then, due to continuity of and t, ^^(t+p; xq, yo) < for 
(xq, yo) close to (xq, ?/o), therefore s{xo, yo) — t{xo,yo) + P, and s is continuous 
at {xo,yo). 

2. Ze{s) > 0, s — Tr. Then z^{t;xo,yo) > for t{xo,yo) + P < t < Tr, thus 
s{xo,yo) = Tr. 

3. Ze{s) — 0. According to Proposition [HI Xs{s) is e-close to x*{t — t) < 0, 
thus Xs{s) < and z^ls) — j{xe{s) + \ze{s)\) < 0. Therefore Zs{t) intersects 
the plane z — transversally. Then, due to continuity of z^, there exists a 
moment s close to s, such that z^{s; xo,yo) — 0. If s < t{xo,yo) + P, then z^ 
decreases for s < t < i{xo,yo) + P, thus Zs{i + p) < 0. 

In all cases, s is continuous. □ 

Proposition 11. Let z^{s£,xo,yo) — with a sufficiently small e. Then both 
v{xo,yo) md u{xs{se),yE{se)) are e-small. 

Proof. According to Proposition[2l is e-small. Then by following the solutions 

We{t) and w*(t) in backward time we obtain that ^^(r) is e-close to w*{t~T), which 
implies that v{xo,yo) is e-small. The second part of this Proposition is obtained 
in the same way by following the solutions We{t) and w*(t) in forward time, as in 
Proposition [51 □ 

Lemma 1. If s is sufficiently small, then operator We{xo,yo) is continuous with 
respect to (xq, yo)- 

Proof. Follows from the definition of Wg, see section[2l and from Proposition[TOl □ 

Lemma 2. Let {x,y) £ Il{a) be a fixed point of the operator We{xo,yo)- Then the 
solution We{t; x, y) is periodic. 

Proof. If suffices to show that if {x, y) is a fixed point, then Case 1 from the defini- 
tion of holds, so that the value of We corresponds to a point on the trajectory 
■We(t;x,y) and is not adjusted, as in Cases 2-4. Note that in Cases 3 and 4, 
Ws{x,y) ^ n, thus they cannot hold. It remains to prove that Case 2 also cannot 
hold. 

Suppose that Case 2 holds, and let, for example, Se{x,y) > a + a. Denote 
{x,y) = We{se,x,y). According to Proposition [Til both v{x,y) and u{x,y) are 
e-small. Moreover, according to Propositions [3 and [SI 

a < Se{x, y)-(T < u{x, y) - v{x, y) + p. 
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Remember that \u{x,y)\ < a/2, and p < a/4. Thus, 

v{x, y) < —a/4. 

Note that the hnear combination in the definition of moves the point {x, y) 
in the direction of the point (0, —2a) in {u, u)-coordinates, thus further decreasing 
v{x,y). Therefore, we obtain v{x,y) < —a/4, which contradicts Proposition 1111 
This proves that only Case 1 can hold for {x,y). □ 

Lemma 3. For sufficiently small e the rotation "f{id — W^, n(a)) of the vector field 
p — We{p) at the boundary of the set n(a) is defined by 

7(id-iy„n(a)) =sgn(A). 

Proof. Let for example, 

A>0. 

Consider and (5_|-, the upper and lower sides of the parallelogram n(a). Then, 
by definition, 

Weix,y) = w*{a + 2a), {x,y) e and 
Weix,y)^w;{a~2a), {x,y) e Q+. 

In other words, 

u{We{x,y))^2a, {x,y)eQ-, (21) 
u{W,{x,y)) ^ -2a, {x,y)eQ+. (22) 
Also, according to Proposition |S1 

\imv{We{x,y))^Q. (23) 

The relationships ([2T |) - ((23l) imply that in the coordinates (u, v) the mapping 
on the boundary of 11 is close to the linear mapping Li{u,v) = (0, — 4w), thus the 
mapping id — is close to (and therefore co-directed with) the linear mapping 
L2{u, v) = {u, 5u), and the result follows from the properties of the rotation number. 

□ 

Lemma [3] implies that the operator W^{xo, yo) has a fixed point (i, y) on the set 
n(a), which defines a periodic solution of ([2]) according to Lemma [2j Moreover, 
according to Propositions [S] and HI v{x, y) and u{x, y) are both e-small, which 
implies both that {x,y) is e-close to {x*,y*), and that the minimal period of the 
solution starting from {x,y) is e-close to cr — t. Thus, Theorem [1] is proved. 

4. Example. Consider the system 

X = —ay + z/a, 
y = X + 1, 
ez — X + \z\. 



In this case (O takes the form 



and (O turns into 



X — —ay + x/a, 
y ^x + 1, 

X = —ay — x/a, 
y — X + 1. 
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For a > oq ~ 1.9 the corresponding curves Wa{t) and Wr{t) intersect transversally 
on the plane {x, y), see Fig. |4l 



y y 




Figure 4. Trajectories Wa{t) (sohd) and Wr{t) (dashed); left: a = 
2, right: a = 3. 

For any a > oq this system has a periodic canard. The last figure graphs the 
numerical approximation of such trajectory, together with the limiting curve, which 
consists of Ta, F^, and a vertical segment connecting them. 




Figure 5. Periodic canard for a = 3 with e = 0.1 (thick) and the 
limiting curve. 
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